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Abstract 

Bardeen's transfer-hamiltonian method is applied to magnetic tunnel junc- 
tions having a general degree of atomic disorder. The results reveal a close 
relationship between magneto-conduction and voltage-driven pseudo-torque, 
and also provide a means of predicting the thickness dependence of tunnel- 
polarization factors. Among the results: 1) The torque generally varies with 
moment direction as sin 9 at constant applied voltage. 2) Whenever polar- 
ization factors are well defined, the voltage-driven torque on each moment is 
uniquely proportional to the polarization factor of the other magnet. 3) At 
finite applied voltage, this relation implies significant voltage-asymmetry in 
the torque. For one sign of voltage the torque remains substantial even if the 
magnetoconductance is greatly diminished. 4) A broadly defined junction 
model, called ideal middle, allows for atomic disorder within the magnets and 
F/I interface regions. In this model, the spin- (a) dependence of a basis- 
state weighting factor proportional to the sum over general state index p of 
(/ / dydzty Pt(7 ) 2 evaluated within the (e.g. vacuum) barrier generalizes the lo- 
cal state density in previous theories of the tunnel-polarization factor. 5) 
For small applied voltage, tunnel-polarization factors remain legitimate up to 
first order in the inverse thickness of the ideal middle. An algebraic formula 
describes the first-order corrections to polarization factors in terms of newly 
defined lateral auto-correllation scales. 

PACS: 85.75.-d 
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1 Introduction 



When first predicted, voltage-driven pseudo-torque in magnetic tunnel junctions 
(MTJs) appeared to be a marginal effect [1]. (Sec. 2 explains our use of the 
prefix pseudo- in the term pseudo-torque.) The lithographic scales and resistances 
available in early experimental MTJs appeared too large to permit anything more 
than a very small torque term in the Landau-Lifshitz equation. Resistive heating of 
the MTJ would have limited its possible consequences to only a small voltage-driven 
decrease of linewidth of narrowly-fo cussed Brillouin scattering. (This prediction was 
never tested.) As a result, one could not yet predict anything as remarkable as the 
now well-established magnetic reversal and high-frequency precession observed when 
the resistive barrier is replaced by a metallic spacer. For recent experimental work 
and earlier references dealing with switching and current-driven oscillations involving 
metallic spacers, see Refs. [2] and [3]. 

But in recent years, experimental activity in tunneling magnetoresistance has 
expanded vastly. It is fueled in great part by the experimental discovery of substan- 
tial tunneling magnetoresistance [4] at room-temperature and the resulting intensive 
exploration of non- volatile magnetic memory reviewed recently [5]. A part of this 
activity is the search for junction compositions and deposition techniques which 
lower the resistance to values more suitable for integrated-circuit application. In- 
deed, there now exist very recent experimental reports of current-driven switching in 
MTJs [6, 7]. This development may make possible two-terminal memory elements 
avoiding resort to three-terminal devices using both a metallic spacer for switching 
and a tunnel barrier for reading [8]. 

According to recent reviews of tunneling magneto-resistance [9, 10, 11, 12], em- 
pirical ferro magnet polarization coefficients Pi [i =L,R refer to left and right magnets 
Fj in Fig. 1(a).] measured with FJS junctions having a superconducting counter 
electrode [13] account well for magneto-resistance in FIF junctions. Let the formula 



for current density at constant applied voltage V define the dimensionless coefficient 
i of magneto-conduction. Here 9 is the angle between the moments. (The — sign 
occurs in Eq. (1) because of the convention in Fig. 1 where particle-number current 
is positive for V > 0.) In this article, the coefficient i is more convenient than the 
experimentally preferred low-voltage tunneling-magnetoresistance ratio 



J(V,9) 



Jo(V)[l +l(V)cos9], with J >0 for V > 



(1) 



TMR = (i? A p - Rp)/Rp = 2l/(1 - l). 



(2) 
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The original equation due to Julliere [14] , expressed in our notation by the formula 

i = PlPr, (3) 

enjoys considerable success in interpreting experiments [9]. We find below that when- 
ever i separates this way into two polarization factors characteristic of the respective 
electrode-and-barrier compositions, pseudo-torque expressions having dimensionless 
coefficients r L and r R [See Eqs. (13), (19), and (20) below.], whose simplicity par- 
allels that of Eqs.(l) and (3), hold also. The presence of the same average current 
density Jq{V) in equations both for magneto-current and torque represents a strong 
connection between these two phenomena. 

After the commonalities in Sees. 2 and 3, these mutual relations (Sees. 4 and 5) 
between magneto-conductance and pseudo-torques constitute the first of two parts 
of the present article. The second part (Sees. 6 and 7) is stimulated by the fact 
that theory does not generally support the separability of spin-channel currents into 
the left- dependent and right-dependent factors needed to justify polarization factors 
in the first place. Previous theories attack the question of polarization coefficients 
within the context of real electron structure by considering the transmission of elec- 
trons initially occupying well-defined crystalline-momentum states [15, 16, 17]. They 
posit either complete absence of disorder or special types of disorder only within the 
barrier to legitimize tunnel-polarization factors. The present approach, detailed 
below, complements those works by excluding disorder only from a subregion of the 
barrier. 

Electron scattering, which causes metallic resistivity, abounds within experimen- 
tal MTJ electrodes. A new feature of the present work is to forego altogether 
crystal-momentum quantization within the electrodes. This feature is particularly 
appropriate to contemporaneous experiments relying for electrodes on evaporated or 
sputtered magnetic elements and alloys having high defect concentration [2, 3, 6, 7]. 
Both alloying and structural defects may cause an electron to scatter many times 
within the electrodes before and after it tunnels across the barrier so that initial and 
final crystal momenta are undefined. 

Our elastic-tunneling theory rests on Bardeen's transfer-hamiltonian method (BTM) 
[18, 19], which is applicable to tunneling transitions between thermal baths of electron 
states without any spatially conserved observables. Bardeen defines two sets of basis 
states — one for the left electrode-and-barrier and one for the right electrode-and- 
barrier. Fermi's "golden rule" for transition rates gives the tunneling current. Thus 
our theory of MTJs has broader application than many others, previously reviewed 
[9, 11, 12], which rely on scattering of Bloch electrons. Although more modern than 
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Figure 1: (a) Scheme of magnetic tunnel junction and key to notations, (b) Equiva- 
lent circuit for spin-channel currents and further key to notations. 



Bardeen's method, they must assume defined initial and final momenta. 

Our model of the junction, called ideal middle, excludes disorder only from a 
central geometric slab of uniform thickness w, which may consist of vacuum or pe- 
riodic crystal lying somewhere within the barrier. We find that exact factoriza- 
tion of channel-to-channel current, which leads to Eqs.(3) and (19) below, occurs 
in the limit w — > oo, just as in the case of complete absence of disorder. Further, 
our parametrization of lateral auto-corellation (See Sec. 7) of the Bardeen basis- 
function sets predicts well-defined tunneling-polarization factors for finite barriers 
to first order in w" 1 , which enhances their legitimacy for interpretation of experi- 
ments involving any degree of disorder. Computations and measurements of the 
new corellation-scale parameters £ CT could shed quantitative light on the genesis of 
polarization factors. 

By way of organization, Section 1 is this Introduction and Section 2 shows how 
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spin-channel tunnel currents generally determine voltage-driven torque. Section 
3 uses the BTM to derive the resulting fully general expressions for the magneto- 
conduction, the torques, and the relevant dimensionless coefficients t, tl, and tr. 
Section 4 shows how tunneling-polarization factors and the resulting simple expres- 
sions for i, tl, and tr arise from a formal separability condition. Section 5 addresses 
the expressions for voltage-unsymmetric torque arising from voltage dependences of 
polarization factors. Section 6 demonstrates the separation condition and derives 
the tunnel-polarization factors which arise in the ideal-middle model at w — > oo. 
Section 7 expands the magnetic tunneling properties for finite w and derives a for- 
mula for the first-order w _1 -dependence of tunnel-polarization factors. Section 8 
summarizes and discusses the results. 

2 First currents, then pseudo-torques 

Whenever two ferromagnets are separated by a nonmagnetic spacer, whether a tun- 
neling barrier or a metal, exchange-generated pseudo-torques acting on the magnetic 
moments are attributable to the flow of spin-polarized current. For a fuller discussion 
of the genesis of pseudo(or effective)-torque from the principle of spin continuity, see 
Appendix B of Ref. [20]. Essential is the interpretation of magnetization dynamics 

(M = dWL/dt) governed by the additive terms in the macroscopic Landau-Lifshitz 

equation. Ordinarily M represents the precession in place of electron-spin mo- 
mentum localized to a volume element dV due to local causes like magnetic field, 
spin-orbit coupling, etc. But the term describing externally driven spin transfer is 
transparently different. It reflects directly the flow of spin momentum directly into 
dV. 

Indeed, the same may be said about the phenomenological exchange stiffness de- 
scribed commonly by the effective field 2AV 2 m, with M =M s m. This truth is masked 
by the derivability of ordinary exchange torque from variation of the stored energy 
density AT li j{dm i / dxj) 2 . Because spin transfer is driven by an externally supplied 
current or voltage, its effect cannot be derived from a stored energy. Therefore, its 
calculation requires direct recourse to spin currents as detailed below. Since this 
distinction between torque and divergence of polarization makes no difference in the 
subsequent application of the Landau-Lifshitz equation, the prefix "pseudo-" will be 
omitted in the remainder of this article. 

Return now to our problem of spin-transfer torque created by external voltage 
applied to the MTJ. Consider particularly the series electric circuit in Fig. 1 (a) in 
which an external voltage V causes electric-current density J to flow in series through 
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a left metallic ferromagnetic film Fl, a thin insulator I serving as a tunnel barrier, 
and finally a grounded right metallic ferromagnetic film F R . By assumption, Fl is 
sufficiently thin for the direction of spontaneous magnetization Ml(i) = — Ml(x)1 
within Fl not to depend on the plane-perpendicular coordinate x; similarly Mr(s) = 
— M R (a;)r within F R . But the spontaneous magnetizations M L and M R may vary with 
x. (Here the three-dimensional unit vectors 1 and r include the angle 9 = cos -1 r • 1.) 
Thus we lay aside those possibilities of forward spin- wave excitation [21] and volume- 
intensive torque [22] , arising from dependence of magnetization direction on x, which 
become significant for larger film thickness and current density. 

One goal is to calculate the component T R of interfacial torque vector T R per 
unit area, acting on M R , which lies orthogonal to r within the instantaneous plane 
common to 1 and r as indicated in Fig. la. (The orientation of the magnetic space 
spanned by 1 or r is completely disconnected from that of position space x, y,z.) A 
general expression for T R [20, 23] reads thus: 

T R = H[J h , + - J L ,_ + ( J R ,_ - J R+ ) cos 6]/2e sin 9 (4) 

Here the left spin-channel electric current densities Jl,± = Jl,±1 flow through plane 
A (See below) in direction x and the right J R ± = J R ±r flow through plane B. The 
factor —H/2e converts any electric channel current to one of spin momentum. A 
similar expression holds for the pseudo-torque Tl on the left moment. The torques 
T R and Tl must generally be included in the dynamic Landau-Lifshitz equations for 
the two magnetic films. 

Although previously applied only to all-metallic multilayers, Eq. (4) may also 
be used when the spacer is an insulator. For its derivation, one posits the non- 
relativistic n-electron hamiltonian including, besides kinetic energy, coulomb terms 
due to external voltage and electron-nuclear and electron-electron interactions. In 
addition, one accepts the microscopically-based approximation, defensible in the case 
of Co, that the transverse (to local M) components of conduction-electron spin po- 
larization created at the two internal I/F interfaces decay to zero well within a 
characteristic distance d± ~ lnm [20], which was estimated explicitly for Co/Cu and 
other interfacial compositions by scattering computations [24]. Moreover, in one 
experiment the threshold current for switching of Co by polarized current flowing 
through a metallic spacer is simply proportional to film thickness down to 1 nm, 
confirming that the transverse polarization inside the ferromagnetic film vanishes 
at this scale [25]. Therefore the currents in the left and right magnets must be 
polarized along instantaneous left (1) and right (r) moment axes at depths greater 
than dj_ from the F/I interfaces. Thus our work excludes thicknesses < 1 nm, which 
require special treatment sensitive to atomic layering [26]. 



6 



In the extensive literature on tunneling magnetoresistance involving Fe, Co, Ni 
and magnetically concentrated alloys of these elements with others of lower atomic 
number, there is little indication of spin relaxation at I/F interfaces. Moreover 
experiments at cryogenic temperatures reveal that the distance Am of spin relaxation 
due to spin-orbit coupling for the polarization component along the axis M is about 
50 nm for Co and about 5.5 nm for Ni-Fe [27]. Thus it follows that, at least in 
the case of Co where X\\ >> d±, the channel currents Jl,± and J Ri ± should be 
evaluated at the planes A and B lying at the distance d± from the respective F/I 
interfaces. For within the space between these planes one may neglect spin-orbit 
effects and embrace the well-known spin-continuity relation which equates the sum 
of equivalent interfacial pseudo-torques with the net inflow of spin current [20, 23], 
having polarization directions 1 on the L side and r on the R side. In the notation 
of Fig. 1 (a), the statistical average of this equality becomes 

T L + T R = A [( j L> _ _ j Li+ )l + (J R+ _ j R! _) r ] . ( 5 ) 

By our assumed neglect of changes in Ml, we write 1 • Tl = 0. Therefore the scalar 
product of Eq. (5) with 1 eliminates T L and gives Eq. (4) for the magnitude T R . A 
similar equation holds for Tl. 

The above argument neglects a decaying and spatially oscillating transverse cur- 
rent, calculated in certain FNF cases to lie between and ~10% of the incident spin 
current (See Fig. 7 of Ref. [24]). It is likely due to specular interference created at 
the perfect interface assumed in the calculation. Studies of FMF exchange coupling 
in vogue 10 years ago suggest that extremals in the Fermi surface determine the 
wavelength and cause the amplitude to decay with distance. The amplitude will be 
decreased by irregularities at real imperfect interfaces. 

Even in the absence of applied electric voltage (V = 0) an additional perpendicular 
component of exchange pseudo-torque T RJ _ = K\ x r = — T L j_ predicted for MTJs 
[1] is generally related to phenomenological coupling energy —K\ • r = — K cos 6. It 
must also be included in the Landau-Lifshitz equation for the dynamics of magnet 
F R However, in that toy rectangular-barrier MTJ model [1], the (uncalculated) de- 
pendence of T R _|_ and T L j_ on applied voltage occurred only in higher order (a V 2 ) 
than the torque given by Eq. (4) (a V). Moreover, its dynamic effect is relatively 
weaker in structures with coincident easy anisotropy axes and low magnetic damping, 
such as the pillars using metallic spacers experimentally favored for efficient current- 
driven switching [2]. Indeed, steady oscillation excited by a steady electric current, 
such as that observed [2, 3], is possible with T RJ _ = 0, but not in the absence of in- 
plane T R . In addition, the BTM used here does not readily provide this out-of-plane 
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torque. For these reasons, we do not attempt to predict the perpendicular torque 
component in this work. 

3 Magneto-conduction and torques 

Equation (4) effectively reduces the interacting-electron problem of voltage-driven 
torque to the customarily independent-electron problem of spin-channel currents. 
One recently reviewed BTM-based theory of collinear MTJ magnetoresistance [11] 
extends naturally to tunneling between spin channels for general 6. For adaptation 
of the BTM [18, 19] to the MTJ of Fig. la, a stationary basis state \p, a) within the 
electron reservoir Fl is assigned orbital index p and majority /minority spin a = ± 
quantized along axis 1. It satisfies (H+eV — e PjCr )\p, a) = 0, and decays exponentially 
within the barrier, considered semi-infinite in width when defining the basis states. 
Here, H = p 2 /2m + E cr |cr)C/ cr (x, y, z)(a\, where the potential U a depends on spin 
within the ferromagnets according to intinerant-electron magnetism theory [28], but 
not within the barrier. Within Fr, a similar state satisfies (if — e qiCr >)\q, a') = with 
quantization axis r. Because the barrier is assumed to dominate all other resistances 
of this circuit, the spin channels are shown in Fig. 1 (b) as shorted in each magnet 
and/or external-contact region by spin lattice relaxation due to spin-orbit coupling. 
One may disregard spin accumulation and the related distinction between electric 
and electrochemical potentials which are important when a non-magnetic metallic 
spacer substitutes for the barrier [29]. U a includes all elastic terms arising from 
atomic disorder due to alloying, defects, interfacial atomic interdiffusion, etc. The 
state indices p, q simply enumerate the exact eigenstates \p, a), \q,cr') of H in the 
Bardeen basis. Each such state incorporates effects of all multiple elastic scatterings 
without limit. 

Employing the spinor transformation connecting quantization axes 1 and r, the 
transfer matrix element takes the form 

(6) 

Direct extension of BTM [30] to our spin-dependent case gives the expression 

-h 2 

lp,a;qA X ) = 1 ^ dydz(lp p>a d x (p q y - (fq^d^p^), (7) 

where the integral is over unit area for coordinate x lying appropriately (see below) 
inside the barrier. The energies e PjCr and e q ^i may differ only infinitesimally from 
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(p,a\H - e\q,a') = 



7p,+;g,- 
— lp,-;q,- 



COS 7j 

• ft 
. sin 77 



7jv 
7p,- 



. sin ■ 

cos 



the Fermi value e = £p- The hamiltonian H, the left (4> p ,a) and right ((p q ,a') orbital 
wave functions, and these matrix elements (7) are real. 

Only the neglect of cross-barrier overlaps (p, a\q, a') allows use of the Fermi golden 
rule of perturbation theory which is strictly valid for an orthonormal basis. Substi- 
tution of the perturbation (6) into this rule is followed by summation over the initial 
states in an infinitessimal energy band of width eV. Thus the partial electric current 
density flowing between channel a in Fl and channel a' in Fr becomes 

— 1ttp 2 V , 

J*,*' = —f—E^iP, *\H- e F \q, a'f (8) 

at T = K. The ' in ^2 pq imposes the conditions Ep < (£ P ,a,£q,a') < £p + eV. 

Notations in the equivalent circuit shown in Fig. 1 (b) make plain the relations 

Jha = Ja,+ + Ju- , Jr<j> = J+,cr' + J-,cr', (c, o' = ±) (9) 

needed in Eq. (4). The right hand sides of these equations are evaluated from Eqs. 
(6-8). 

Next we write the total electric current density J = Jl,+ + Jl,~- With the notation 

2rreV i 

= —f—Y, P , q il*- q ,„> (io) 

for interchannel particle-number tunneling conduction with the angular factor omit- 
ted, the above equations combine to give Eq. (1) with 

j = e(r +i+ + r_,_ + r +> _ + r_ >+ )/2 (11) 

and the electric magneto-conduction coefficient 

i = e(r +i+ + r_,_ - r + ,_ - r_, + )/2j . (12) 

Eq.(4) becomes 

T R = -(/ir R J /2e)sin^ (13) 

or, in coordinate-free form 

T R = (hr R J /2e)r x (1 x r), (14) 

with the torque coefficient 

r R = e (r +i+ + r + ,_ - r_,_ - r_, + )/2j . (15) 

The fact that the linear combination of the parameters T a y appearing in Eq. (12) 
differs from that in Eq. (15) and a similar one for Tl precludes any fully general 
connection between torques and electrical current. 
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4 Left-right separability and polarization factors 



Particularly interesting relations arise if the summation in Eq. (10) for the inter- 
channel particle current happens to separate into left- and right-dependent factors 
in the form 

IV = f^L,A,a'- (16) 

Here the coefficient / , which we make no attempt to evaluate, is independent of 
a, a'. (Sections 6 and 7 address conditions for this separability.) For then Eq. (11) 
gives 

J = ^(fi L)+ + JV)( fi R,+ + fi R,-) (17) 
and Eq. (12) gives Eq. (3) with the tunneling polarization parameters 

(i = L,R) (18) 



which are directly measurable using FIS junctions [9]. In these terms, Eqs. (1) and 
(3) give the magneto-conduction and Eq. (14) the torque with 

r R = P L . (19) 

Similarly, the torque on the left magnet is 

T L = -(hr L J /2e)sm9, r L = P R (20) 

or, in coordinate-free form 

T L = ^J 1 x (r x 1). (21) 

The Eqs. (3), (19), and (21) show the very close relation between current-driven 
torques and magneto-conduction at the same voltage, summarized by i = TlTr, if 
the separability condition (16) is satisfied. 

The ground-breaking paper of Julliere [14] gave equations equivalent to (3) and 
(18) taking fl^ and Qr ct / to be spin-dependent basis-state densities at e — £p- It 
appeared to attribute the dimensionless magneto-current coefficient i = PlPr to bulk 
properties of the two magnetic compositions involved. But the analytically solved 
free-electron rectangular-potential model [1] shows that an interface-dependent factor 
must be included in as well. The transfer-hamiltonian treatment of this toy 
model follows immediately from the spinless treatment [30] giving 

*V = *W(»o + k t) ( 22 ) 
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where 



k 2 ^ = 2mE it<T /h 2 and Kq 



l = 2mB/H 2 . 



(23) 



Here, E ii<T is the kinetic energy at the Fermi level and B is the barrier potential 
measured from the Fermi level. Equation (18) now gives 



in agreement with Ref. [1]. In this formula, the first factor depends purely on basis- 
state densities in the magnet, while the second mixes magnet and barrier properties. 
The results of the toy model [1] satisfy the general magneto-conduction relations (1), 
(3) and torque relations (14), (19), (21) with this substitution. 

We note in passing that experimental variation of barrier height B shows consid- 
erable support for the zero of i at the barrier potential satisfying — k it+ k it _ = 
expected from Eq. (24) [31] (for small V). Therefore, in spite of its fundamental 
naivete, this toy model enjoys some degree of credibility. It illustrates the general 
fact that, even when separability holds, each polarization factor is a property of the 
electron structure of the magnet and barrier combination as demonstrated by many 
experiments and calculations. Section 7 will discuss how tunnel polarization may 
vary with barrier thickness. 

5 Finite bias and torque asymmetry 

In experiments, TMR typically decreases significantly with increasing finite V [9]. 
Voltage-dependence of interfacial transmission, special state density distributions, 
extrinsic impurity effects, and inelastic tunneling contribute to this decrease [9, 12]. 
This is important because large voltages will be required to read and write in a 
2-terminal memory element. 

The toy polarizations of Eq. (24) will serve to illustrate qualitatively the very 
unsymmetric effect of finite V on voltage-driven pseudo-torque. One calculation of 
TMR uses the WKB approximation for the free-electron wave function within the 
constant-slope barrier potential sketched in Fig. 2 [32]. The interfacial transmissions 
are approximated by those of the flat-potential polarizations (24). The authors cite 
some experimental support for their results. 

It is the decrease of Pi in the particular electrode which collects the tunneled 
electrons that primarily accounts for the decrease of i in the calculated result [32]. 
In Fig. 2, for V > 0, the collecting electrode lies on the right. Note that the electrons 
whose energy lies in a narrow band (shaded in Fig. 2) just below the Fermi level 




(24) 
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of the emitting electrode on the left of the barrier dominate the tunneling current 
because of the strong energy dependence of the WKB factor exp[— 2 J K(x)dx\ in 
the transmission coefficient. Since these hot electrons lie an amount well above 
the Fermi level on the right, this energy shift eV must be taken into account when 
estimating P R . 




Figure 2: Schematic junction potential for finite V. The shaded bar indicates the 
energy range of most of the tunneling electrons. 



We simplify this model one step further and neglect the width of the shaded 
current band in Fig. 2. It is then clear that Eqs. (23) and (24) with i =L are 
still correct for P L , neglecting correction for the finite slope of the barrier potential. 
However, the equations 

k 2 Ra = 2m(E Ra + eV)/h 2 and k\ = 2m(B - eV)/H 2 , (25) 

obtained by adding eV to each electron energy on the right, must replace Eqs. (23) 
for i =R. 

Figure 3 plots the curves tl = Pr and tr = Pl evaluated from the preceding three 
equations as well as TMR from Eqs. (2) and (3) versus V for the special example of 
a symmetric junction with the parameters /cl- = &r- = k_, = k R+ = 10/c_, and 
k = 6.4/c_ i whereby each electrode has the V = polarization Pl=Pr = 0.5. In this 
illustration, TMR(V) is symmetric because it involves both P L and Pr but Pl,r(V) 
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and the torque coefficients tl,r(V) are not. Although the theory preceding this 
section assumed small V, the present discussion makes reasonable the application of 
the results to finite V with the understanding that the polarization of the collecting 
electrode generally depends more strongly on V. Of course, this toy calculation 
cannot make quantitative predictions of the ^-dependence which must rest on details 
of electron structure [9, 12]. 



0.7 i 




-4 -2 2 4 



voltage V (arbitrary units) 



Figure 3: Schematic effect of finite voltage on TMR, polarization, and torque coeffi- 
cients illustrated by the toy free electron model of a physically symmetric magnetic 
tunnel junction. Note that TMR is symmetric but the other coefficients are not. 
The parameters are k q = 6. 4k-, k + = 10k-. 



Note that while critical current density for magnetic excitation is appropriate to 
junctions with metallic spacers, the high resistance of a MTJ makes critical voltage 
more appropriate. (Indeed, strictly speaking, the critical current of a constant- 
current generator will generally differ from the current density flowing at threshold 
in the presence of constant external voltage.) Another significant difference between 
metallic and insulating spacers lies in the angular symmetry of the torque. The 
fixed sin ^-dependence at constant V in the tunneling case has no counterpart in the 
metallic case where more general torque expressions typically contribute to asym- 
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metry of excitation threshold [20]. Now we see that the non-ohmic resistance of a 
tunneling barrier gives rise to the torque asymmetry of tr(V) exhibited in Fig. 3, 
which naturally reflects in yet another origin for asymmetry of voltage threshold. 



6 Ideal- middle model for separability 



A recent publication compares existing theoretical arguments supporting the exis- 
tence of tunnel-polarization factors [17]. Each of them assumes incident states with 
definite crystalline momentum. One common type of argument assumes complete 
absence of disorder so that the tunneling through a thick barrier is dominated by a 
single value of lateral momentum. A different model of Tsymbal and Pettifor [15] 
recovers factorization and therefore the Julliere formula in a tight-binding single- 
band model disordered only within the barrier. Similarly, the model of Mathon and 
Umerski attributes the factorization to phase decoherence due to disorder within the 
barrier [9, 16]. These treatments are augmented with arguments based on the Feyn- 
man path integral in a disordered barrier [17]. Our treatment below complements 
these arguments with the contrary tack of foregoing lateral momentum quantization 
completely within the electrodes and I/F interfaces while preserving ideal crystalline 
ordering or vacuum within the middle of the barrier. 

Figure 4 indicates the structural scheme. The left (ip p>eT ) and right (<p g , a ') orbital 
basis functions for the transfer matrix, introduced in Sec. 3, are governed in detail by 
the general potential U a or a > depending on crystal structure, alloy composition, de- 
fects, F/I interface roughness and atomic interdiffusion, etc. The quantum numbers 
p and q do not refer to any diagonal operator. Exceptionally, the ideal-middle B of 
the barrier consists of an ideal crystalline slab or vacuum region defined by a < x < b 
where the planes x—a, b are dubbed portals of the ideal middle. In order to define 
the left and right basis-state sets of the Bardeen theory, the barrier potential extends 
into alternative semi-infinite spaces (a < x) and (x < b) , where it is greater than 
£f, independent of or periodically dependent on y, z and independent of a and a'. 
The respective conditions ip PtU — > for x — > oo and ip q y — > for x — > — oo complete 
the definitions of ip p ^ and <p q y. 

The effective-mass theorem [33] is valid when e is near the bottom k = k of the 
conduction band within region B. Then the evanescent portion of a left-magnet 
basis function within this region is approximated by 



where ^ P}C satisfies (7ibar — £ P ,a)^ P ,a=0 an d ^ p>a — > for x — > oo, and Wcb.ko is 
the Bloch function for the bottom of the conduction band. The effective barrier 




(26) 



14 



w 




a b 



Figure 4: Depiction of the ideal-middle model of a magnetic tunneling junction. 
Disorder without limit is permitted in both electrodes and barrier except within a 
central slab B of the barrier lying between the portal planes x — a,b. 



hamiltonian is 7Ybar = —h 2 V 2 /2m c ^ +U(x) where m c b is the effective mass and 
U(x) (> £f) is the spin-independent atomically smoothed effective barrier potential. 
Similarly for Fr, <p q y = ^q^'Ucb,^ with § q y —> for x — > — oo. In case of vacuum, 
(\l/, $) are indistinguishable from (ijj,(p). (Note however that this treatment fails if 
both V is finite and the FI interfaces are disordered, for then U depends on y and z 
as well as x.) 

Assuming periodic boundary conditions in the s= (y, z) sub-space, the evanescent 
portions of left and right basis states within B are conveniently fourier-expanded in 
space s with the WKB approximation giving 



and 



r px 

~ \>,o-(k)[ft(&, a)/n(k, x)} 1 ^ 2 exp — / n(k, x')dx + ik • s 
k I Ja 

®g,a' — / i g,cr'(k)[/«(^, b)/tz(k, x)} 1 ^ 2 exp — / K,(k, x')dx' + ik • s 

k . J x 



(27) 



(28) 



where the sums are carried over a 2-dimensional reduced Brillouin zone. These 
formulas employ the function 
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K,(k,x) = [kI(x) + k 2 ] , with k\ = 2m ch [U (x) - e F ]/h 2 , (29) 

where m is the imaginary component of the wave-vector in region B. Note that Eqs. 
(27) and (28) reduce to expansions of ip P:(T and </? ?;(T / with coefficients A PjCr (k) and 
/i g)(T /(k) on the portal planes x = a and x = b respectively. 

The transfer-hamiltonian matrix element of Eq. (7) is evaluated at any x lying 
within the interval a < x < b. Consequently $ , and m c b may replace ip, ip, 
and m respectively in this formula. One convenient choice to evaluate Eq. (7) 
is x = x max , satisfying U(x) < U(x max ) for all x, because the resulting condition 
dK /dx(x max ) = simplifies the mathematics. (Inclusion in U of the image potential 
due to electron-electron correlation will often insure the presence of a maximum, even 
if \ V\ is large.) Substitution of Eqs. (27) and(28) followed by integration over y and 
z, with the assistance of the identity J ds 2 exp[i(k — k')-s] = <5 kk ,, reduces Eq. (7) to 

lp,*-, q y = S k F(w, k)A* )<y (k)^, iff /(k) (30) 

where 

F(w,k)= -/t 1/2 (k,a)K 1/2 (k,6)exp[- / dxK(k,x)]. (31) 

Here we use the barrier-middle thickness w = b — a, and note A* CT (k) = A P)<T (— k) 
and //* ,(k) = i^ q y(— k) because ^/ PjCr and are real. [When iu varies in our 
discussion below, A P)(T (k) and /x 9)0 ./(k) remain unchanged because they pertain to 
the semi-infinite barrier independent of w. We merely expand or contract the ideal 
middle of the barrier in Eq. (31).] After rearranging the order of sums, Eq. (10) 
with substitution of (30) and (31) becomes 

IV — ^k^^' k) 5^ k , F{ w i k')Ar(k, \s.')M a '(k; k') (32) 

where each of the two functions 

C a = ^,A^(k) VW, = E^(k)^(k') (33) 

depends only on parameters of the left and right magnet-and-barrier combinations 
respectively. The ' on S' signifies the conditions given previously for Eq. (8). 

In the presence of atomic disorder, the sums in Eqs. (33) are carried over many 
states of randomized character. Therefore they have the nature of statistical auto- 
correlations in y, z-space which should depend smoothly on k and k' and are Taylor- 
expandable about k = k' = 0. (See Sec. 5 for the very different toy free-electron case 
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of vanishing disorder [1], in which one may formally replace p — >k", q ^k w so that 
C a and Ai a ' become proportional to S k y.) In addition, with increasing thickness 
w = b — a of region B, the exponential in Eq.(31) becomes ever more sharply peaked 
at k = 0. Summation over k and k' of the terms in these Taylor series' for finite w 
gives the corresponding terms 

r^H = CH + r«H+ ... (34) 

The initial constants in both Taylor expansions yield 

CW = /W<i< (35) 

with = £ CT (0,0) and fi^, = M. a '(0, 0). Here factors independent of a and a' 
are absorbed into /. Therefore, to leading order in this expansion, the integrations 
in Eq. (32) tend to the left-right separation of the form (16). 

Written in full, the parameters needed in the general polarization formula (18) 
are, to lowest order in the Taylor expansions of Eqs. (33), the basis-state weights 

Q S = E (/ / dydz^{a, y, z)) (36) 

Q t' = E (/ / dydz* 9 Ab, V, *)) (37) 

where J J dydz is carried over unit junction area at the portal positions a and b. 

[See the next section for development of r^,(w).] Note that the latter two equations 
differ generally from the local state (or charge) density often cited in connection with 
tunneling. (LSD oc J J dydz^ ) They reduce to the LSD in the complete absence 
of disorder when each of the two sums reduces to a single term ^^ =Q a and $k=o,o-' 
independent of y and z. 



7 Correction of polarization at finite thickness 

The non-orthogonality between left and right basis functions constitutes a basic 
weakness of the BTM. Even though the validity of golden-rule transition rates 
in BTM is not generally assured, it has an enormous acceptance in the literature. 
The toy free-electron MTJ theory, though founded directly on a solution of the wave 
equation in the entire ideal non-disordered FIF system having a flat barrier potential, 
was evaluated only to leading order in the exponential parameter e~ KW [1]. The BTM 
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calculation for the same model agrees exactly with its results, as one knew it should 
from previous non-spin dependent tunneling theory [19] . 

Let us assume that BTM is correct to the same exponential degree for our ideal- 
middle model as for the toy model. The previous section showed that the BTM 
supports tunnel-polarization phenomenology in lowest order. Continuing with BTM, 
we derive here a correction to the formulas (18), (36), and (37) for polarization which 
we find below varies algebraically, not exponentially, with w" 1 . Therefore these 
corrections should be reliable in spite of this general weakness of the BTM. 

Further progress requires parametrization of the autocorellation functions defined 
by Eq. (33). Note first the consequence of assuming that the possibly disordered 
atomic configuration in F L produces no electrostatic potential in F R and vice versa. 
From Eqs. (27), (28), and (33), in-plane translation of the (disordered) microscopic 
potential of only the left electrode according to s — > s + (B,C), where (B,C) is 
a periodic-lattice translation of the barrier middle, has the effects, from Eq. (33), 
L a — > L a exp[i(k'— k) ■ (B,C)] and M. a > — > M. a '- Averaging over all possible such 
phase changes makes C a and M a > diagonal and eliminates all terms with k ^ k' from 
the double sum in Eq. (32). This equation now becomes 

IV = ^£ k F 2 (k)£ CT (k)A^k) (38) 

using the now diagonal forms of C a and Ai a >. 

Parenthetically, note that in the special case of vanishing disorder, the state 
indices p and q become m, k and n, k respectively, with m, n the respective band 
indices and k the lateral crystalline momentum. Let the basis states be normalized 
to unity. Then the diagonal elements of formulas (33) reduce to 

£ CT = E m |AV k )| 2 M,V k ) , M a > = En, CT '|/W'( k )lV^V k ) ( 39 ) 

with factors independent of a and a' omitted. Here v XymyCr = de mia (\t)/dk x and 
Vx,n,a' = de n ^'(k)/dk x are velocity components normal to the junction plane. Their 
presence in these formulas follows from the restriction on E' in the basic formula (8). 

To evaluate Eq. (38) for finite disorder, specialize to small V and constant U 
inside B. After evaluation of the integral in Eq. (31), it reduces to the form 

IV = /i^ k ^(k)e- 2 - K ( k )£ CT (k)^ CT ,(k) 

where /i does not depend on a or a'. For large w, this sum weights small k heavily, 
as mentioned above. Therefore parametrize C a and Ai a ' for small k with the lateral 
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spatial correlation scales (£<t,V(t') defined by the formulas 



£ a (k) = £ a (0)[l - gk 2 + <D(k% (40) 
M a ,(k) = M a ,(0)[l - r)l,k 2 + 0{k% 

and approximate Eq. (29) with k ~ + (k 2 /2kq) in the exponent of Eq. (31). After 
approximating (over one BZ) with an infinite integral, one finds by elementary 
integration a result equivalent, to first order in w^ 1 , to 

vh^Ch+Ch ( 41 ) 

- / 2 M£ CT (0) (l - -^4t) A*„(0) (l - -^4t) (42) 

where once again factors independent of both a and a' are absorbed into / 2 . Thus to 
this approximation, once again has the factored form (16). (It appears that in 
order w~ 2 , does not separate this way into left- and right-dependent factors.) 
The corrected left polarization factor, according to Eq.(18) reduces upon expansion 
to 

P L = pW + l(l- 2 ) Koie ~ ^ } + - with (43) 
2 V / w + k q 

(o) £ + (0)-£-(0) 

L " £+(0)+£_(0)' 

and similarly for P R . Thus, from given Bardeen basis functions one can obtain polar- 
ization factors, correctly to order w -1 , in a disordered electrode-barrier combination. 



8 Discussion 

Although it is valid only in the limit of weak transmission, predictions from Bardeen's 
tunneling theory [18] are interesting because it does not require electron momentum 
within the electrodes to be conserved. Our application to elastic tunneling through 
ordered or disordered magnetic tunneling junctions yields these new conclusions: 

• In Section 3 we found that the torque at constant external voltage is gener- 
ally proportional to sin# [Eq.(13)] . This result is a direct consequence of the 
single-transition nature of tunneling and the simple form of the spinor transfor- 
mation (6). It contrasts with the more general angular dependence conditioned 
on electron structure and spin-channel resistance parameters in the case of a 
metallic spacer [20]. 
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In general, polarization factors do not exist in the absence of special assump- 
tions, in agreement with previous theory [9, 12]. 

In Section 4 we found that if the polarization factors are well defined, then 
at constant applied voltage the electric current and in-plane torque obey the 
relations (1), (13), and (20). These similar relations are inter-connected by the 
presence of the common factor Jq(V) which we do not attempt to calculate. 
The dimensionless coefficients in these relations are expressed in terms of the 
polarizations by r R = P L , r L = P R , i = PlPr, implying i = r L r R . In particular, 
these general relations are satisfied by the special results of a direct solution 
of the Schroedinger equation for the toy model of parabolic bands and ideal 
rectangular potential barrier [1]. 

Experimentally, TMR is known to usually diminish with increasing external 
voltage V [9, 10]. In Section 5 we considered that it is the polarizing factor 
of the collector electrode which decreases more strongly with V, resulting in 
the unsymmetric schematic pattern of voltage dependence of torque indicated 
in Fig. 3. This lack of symmetry due to the relations tr = Pl and tl = 
Pr implies that the threshold voltage for initiation of dynamic excitation will 
be increasingly asymmetric at the higher values (>100 mV) likely needed for 
writing in memory. Cases may well arise in which voltage-driven switching 
works in only one direction. For selected experimental junctions, switching is 
observed at a voltage sufficiently high for TMR to become negligible [7]. Our 
Fig. 3 indicates how this may happen for switching in but one direction, from 
AP to P. However, our theory would not explain any symmetric persistence 
of switching at voltages great enough to destroy TMR, if this is observed. 

Our approach to the validation of polarization factors complements previous 
studies which accounted for atomic disorder in the barrier assuming electrode 
states with well-defined crystalline momentum [15, 16, 17]. We assume that 
the barrier is thick and includes an ideal crystalline or vacuum middle region 
of thickness w as in Fig. 4. Then a newly derived polarization factor, given 
by Eq. (43), is valid to first order in w~ l even in the presence of disorder 
in the electrodes and interfaces sufficient to destroy the conservation of lateral 
crystalline momentum throughout the electrode and interface regions. The key 
basis-state weight factors (36) and (37) are more general than the conventional 
local state density. 

Our conclusion that the validity of polarization factors increases with increasing 
w tends to undermine our predictions of voltage asymmetry of torque shown 
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schematically in Fig. 3. For, experimental spin-transfer effects such as switch- 
ing will require very thin barriers, making the separability condition assumed 
in Fig. 3 less valid. Previous proposals [12, 17] that validity of polarization 
factors is attributable to unique defect states or amorphicity in the barrier are 
more promising in this respect. 

• Belashenko and co-authors [17] find that certain first-principle TMR compu- 
tations for realistic barrier thickness may be poorly approximated by propor- 
tionality to e~ KW . This casts additional doubt on the applicability of the ideal 
middle to the very thin junctions needed for spin-momentum transfer experi- 
ments. However, our conclusions from this model may bear significantly on 
magneto-resistance experiments carried out with greater thickness, as suggested 
below. 

• Our parametrized expression (43) for dependence of tunnel polarization on 
ideal- middle thickness w is without precedent. A strong dependence is ex- 
pected from certain compositions, like Co, Ni, and certain alloys, such as 
FeCo, lying on the negative-slope region of the Slater-Pauling curve [34]; for, 
their strong contrast between heavily 4sp- weighted density of majority-spin and 
heavily 3d-weighted density of minority-spin bands may be reflected in strongly 
contrasting magnitudes of left lateral autocorellation scales £ + and Theo- 
retical estimation of the left polarization factor will require prior first-principle 
computation of the Bardeen basis functions ip PtCT for the disordered electrode- 
barrier system. From these, one must invert the series (27) to evaluate the 
diagonal elements of the Fourier coefficients X Pj(T . Then Taylor expansion of 
the diagonal element in the first Eq. (33) for substitution into the first Eq. 
(40) provides the coefficients £±(0) and £±. These parameters must then be 
substituted into Eqs. (43) to obtain the left polarization factor. 

• In fact, experimental junctions having composition Fe/Al 2 03/FeCo show de- 
pendence of TMR on barrier thickness [10] at T=2 K where our assumption 
of elastic tunneling should be valid. A monotonic dependence on thickness, 
expected from Eq. (42), is observed for two crystallographic orientations on 
single-crystal Fe, but not on the third. Although the (say) right electrode 
(FeCo) lies on the negative-slope side, the left electrode (Fe) lies on the positive- 
slope side of the Slater-Pauling curve where high 3d density exists for both signs 
of spin so that there may be little difference between £ + and Junctions 
with both electrodes taken from the negative-slope side may yield a more pro- 
nounced thickness dependence of TMR on barrier thickness according to the 
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present theory. 
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